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The bilayer systems exhibit the Bose-Einstein condensation of excitons that emerge due to 
D , Coulomb pairing of electrons belonging to one layer with the holes belonging to the other layer. 

^"*^ ■ Here we present the microscopic derivation of the dynamic equation for the condensate wave func- 

tion at a low density of electron-hole (e — h) pairs in a strong magnetic field perpendicular to the 
layers and an electric field directed along the layers. From this equation we obtain the dispersion 
law for collective excitations of the condensate and calculate the electric charge of the vortex in the 
exciton condensate. The critical interlayer spacing, the excess of which leads to a collapse of the 
superfluid state, is estimated. In bilayer systems with curved conducting layers, the effective mass 
of the e — h pair becomes the function of the e — h pair coordinates, the regions arise, where the 
' energy of the e — h pair is lowered (exciton traps), and lastly e — h pairs can gain the polarization 

in the basal plane. This polarization leads to the appearance of quantized vortices even at zero 
temperature. 
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I. INTRODUCTION 
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In recent years considerable attention has been focused on the so-called bilayers, i.e. the systems that represent 
two thin and closely-spaced conducting layers separated by the energy barrier for carriers^. This has been mainly 
O \ caused by a successful realization of such systems in GaAs/AlGaAs heterostructures. At small interlayer spacings 
the Coulomb interaction of carriers belonging to adjacent layers proves to be rather essential. The interaction as well 
as a two-dimensional character of carrier motion determine the unique properties of bilayers. In particular, in the 
bilayer system there may arise a collective state, which is characterized by the interlayer phase coherence and, as a 
consequence, the system becomes superfluid. 

One of the examples of the bilayer system is the n — p system consisting of the n-layer with electron conductivity 
and the p-layer with hole conductivity. In the development of the ideas by Keldysh and collaborators^ it has been 
established that in this system in the weak-coupling limit the attraction of electrons from the n-layer conduction band 
and holes from the p-layer valence band gives rise to the coherent BCS-type state of electron-hole (e — ft) pairs^. By 
\ virtue of the spatial separation of pair components, the superfluid motion of e — ft pairs reveals itself in an emergency 
of equal and opposite directed dissipationless electric currents, those may be detected in the experiment. 

The pairing of electrons and holes as well as the condensation of e — ft. pairs is favored in a strong magnetic field 
H, normal to the layers, provided that the magnetic length Ih (Ih — \J hc/eH) is considerably shorter than the Bohr 
radii of both the electron and the hole, and the carriers in the both layers occupy the states in the lowest Landau 
level^. In this case (in the symmetric gauge), the wave functions of both the electron and the hole are "localized" 
in the region of about Ih in size (rather than the size of the effective Bohr radius), that substantially increases the 
\ binding energy of the e — ft pair at small interlayer spacings (d < Ih)- Note that at low pair densities, n p <C l H 2 , the 
• • ■ electron-hole pairing occurs in the real space (strong-coupling limit of e — ft pairs) rather than in the momentum space. 
The properties of an individual e — ft pair in the two-layer system (indirect magnetic exciton) have been considered 
in Ref.6. 

Another example of the phase-coherent system is the n — n system including two electron layers in the magnetic 
field perpendicular to the layers with the total filling factor Vt = V\ + v% = 1. Despite the apparent difference between 
the n — n and n — p systems, they are, in essence, fully similar. Indeed, at the equality i>t = 1, the number of electrons 
in the first layer is equal to the number of free states in the lowest Landau level in the second layer (for certainty, 
we assume v\ < 1/2). Then the transition to the "hole" representation in the second layer— gives the one-to-one 
correspondence between the n — n system and the above-discussed n—p system in the magnetic field. So, in the 
systems of both types the interlayer phase coherence is the consequence of Bose condensation of e — ft pairs, and due 
to this reason the superfluidity in bilayer systems is often called the " exciton superfluidity" . 

As mentioned above, the exciton superfluidity phenomenon was predicted for the n—p systems in the weak-coupling 
limit and at a high density of e — ft pairs rather long ago^. At the present time, one has impressive evidence that 
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the superfluidity in question was revealed experimentally in quantum Hall n — n systems at v\ — ^ 2 — 1 /^i 9 -. This 
experimental situation corresponds to the intermediate pair density, when the distance between the pairs is of order 
of the pair size. We stress that the v\ = v% = 1/2 condition is not necessary: for the exciton superfluidity the equality 
between electron densities and hole densities, i.e., Vt = 1, is necessary. 

The experimental revealing of exciton superfluidity in bilayer systems is a rather strong incentive to its further 
studies. The study of superfluidity of strongly bound e — h pairs in the low-density limit is of particular interest from 
the viewpoint of theory, because one can expect that the e — h pairs with spatially separated components behave as 
a weakly nonideal Bose gasiS, which is described by a relatively simple differential equation of the Gross-Pitaevskii 
typeii. In the present paper in the mean field approximation we derive microscopically this equation for the condensate 
wave function. We consider the low-density limit of e — h pairs in quantum Hall n — n systems with — 1. It is 
assumed that the magnetic length is considerably shorter than the Bohr radii of electrons and holes. Besides, it is 
considered that the interlayer tunneling is negligibly small. The equation is applied to obtain the dispersion law of 
collective excitations and electric charge of the vortex in the condensate of electrically neutral e — h pairs. 



II. THE DYNAMIC EQUATION 



We now proceed to the derivation of the dynamic equation. The Hamiltonian of the bilayer electron-hole system 
in the perpendicular magnetic field is 



H= ( dvY J i>t(?)hk(v)M*) + \ f rfrdr'^^(r-r')V'fc(r)V i + (r')^(r')^(r) 

J k k,l 



(1) 



Here and further on it is assumed that the indices fc, I = 1(2) indicate the electron (hole) layer. The operators ^(r) 
and V'/c( r ) ar e the operators of electron (hole) creation and annihilation at the point with two-dimensional radius- 
vector r. The spin indices are omitted, since electrons are considered spin-polarized. The kinetic energy operator has 
the form 



hk(r) 



2m k 



or c i 



(2) 



and the Coulomb interaction energy of carriers is given by 



V kl (r-r') = {-l) k+l - e " ==. (3) 

eV(r-r') 2 + (l-^)d 2 

In formulas is the electron mass (fc = 1) or the hole mass (fc = 2), e is the dielectric constant, d is the 

interlayer spacing, 5k,i is the Kronecker delta-symbol. 

To derive the dynamic equation for the condensate wave function we shall follow Keldysh's paper—, according to 
which the coherent state of e—h pairs, in the mean field approximation, may be represented by the vector \<f)) ~ D l p\Q) 1 
where the unitary operator is written as 



D,^ = exp{ J dr 1 dr 2 [^(r 1 )$(r 1 ,r 2 ,i)e^ t ^+(r 2 )- V 2 (r 2 )$*(ri,r 2 ,i)e v Vi(ri)]}. (4) 

The vector |0) is the vacuum state of the system: ^|0) = 0. The unknown function <E>(ri,r 2 ,£) and the chemical 
potential fj, of e — h pairs, which enter into D^, should be found from the Schrodinger equation for the vector \<f>) . 
This equation is conveniently written as 



(ihD+^D^ - D+HD^O) = 0. (5) 

In the expression D^HD^ the operators D~t and D$ realize the linear transformation of the creation/annihilation 
operators of electrons and holes. For example, the operator ipi (r) is transformed as 
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where with an accuracy up to terms cubed in <I> we have 



C(r,r') = J(r-r') - (1/2) / dr"$(r, r")$*(r', r"); 



5(r,r / ) = $(r,r')-(l/6) / / dr"dr'"$(r, r")$*(r"', r")$(r"', r'). 
As a result ol transformations, Eg .([5]) in the Hartree-Fock approximation takes the following form: 

f drxdr^^rOMr!,^,^ (6) 
J k 

We do not give here the explicit expressions for h k (ri,r 2 ,t) and <2(ri,r 2 ,t) because they are too cumbersome (for 
more details, see Ref. 12). Since -0fc|O) = 0, equality flBJ is equivalent to the equality Q(r 1 , r 2 , t) = 0, that represents a 
nonlinear integro-differential equation for the function $(ri,r 2 ,i). In the general case, Q comprises the convolutions 
over the coordinates of an arbitrary large odd number of functions However, in what follows we shall retain the 
terms not higher than the third order in $, because we are interested in the low-density limit of e — h pairs. We first 
consider the non-interacting pairs. For this purpose, in the equation Q(ri,r 2 ,i) = we keep only the linear in <f> 
terms. We obtain 



ih— = [Jn(n) + h 2 (r 2 ) + y 12 (n - r 2 )]$ - n§. (7) 

Equation (J7]) can be transformed into the equation for the wave function of the e — h pair condensate (slowly varying 
in space and time) by averaging it over "fast" variables that describe the intrinsic degrees of freedom of e — h pairs. 
To separate the intrinsic degrees of freedom, we turn from the electron/hole coordinates ri and r 2 to the difference 
coordinate f = ri — r 2 and to the coordinate of the center of inertia of the pair R = (rriiYi + m 2 r 2 )/(mi + m 2 ). In 
the symmetric gauge A/.(r/.) = (1/2)H x , it is convenient to seek the approximate solution of Eq.([7]) for the states 
belonging to the lowest branch of the exciton spectrum in the following form: 



$(f , R, t) = exp(^R • (H x f ) + If • A) ^ o(i . _ ( R) <). ( 8 ) 

Here 7 = (m 2 — mi)/(m 2 + mi), fo{ r ) = {V^Ih ) 1 exp(— r 2 /4Z^) is the wave function of the state at the lowest 
Landau level, and p = —ij^(H X ^-). The structure of expression © is determined by the solution of the steady- 
state Schrodinger equation for the e — h pair in the bilayer system. Indeed, at \&(R) oc exp(^p • R) expression (j8]) 
turns into the wave function of the exciton with the momentum equals to p^^. In the case of a slowly varying 
function ^(R) expression (JS| presents a wave packet composed of low- momenta exciton states. To obtain the sought- 

for equation for ^(R, t), we act upon Eq.([7]) by the operator (fo(r — p) exp^— jf^R • (H x f) — ■ and integrate 

over the variable f related to the intrinsic degrees of freedom of the e — h pair. As a result of integration, the first two 
terms in the right-hand side of Eq.0 give (?io; c /2)W( R , t), where lo c = (— + r^-) — ■ Since R is the "slow" variable, 
the integral of the Coulomb term in Eq.([7]) can be calculated by expanding ipo(v — p) in powers p to an accuracy 
of the second-order terms. Putting the chemical potential equal to the ground-state energy of noninteracting pairs, 
p = (Hlo c /2) + Eq, where the binding energy of the e — h pair is given by 
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Eq.([7]) can be transformed into the Schrodinger equation for the wave function of the condensate of noninteracting 
e — h pairs: 



_dtf(R,t) h 2 d 2 Tfn . 

ih — 7- — - — - = ■— - *(R,t). (10) 

dt 2M H dK 2 y ' ' v ' 

In the case considered (magnetic length Ih is considerably shorter than the Bohr radii of both the electron and the 
hole), the effective mass Mh of the e — h pair arises due to the Coulomb interaction between the electron and the 
hole, and is independent of their masses m\ and m-Jl^: 




The mass Mh increases with both the increasing magnetic field and the increasing spacing between the layers. 

To take into account the interaction between e — h pairs, it is necessary to include the cubic terms in the equation 
Q(ri, Y2,t) = apart from the linear in $ terms. At a low density of e — h pairs the contribution from these terms 
can also be calculated with the use of the perturbation theory. The correction to the chemical potential of the e — h 
pair is equal to 



(i) 



j[4nd- (2nf/ 2 l H + (2^ 2 l H exp (^)erfc(-^-) 



(12) 



In Eq. (fT2")) the first term is related to the energy of the electric field between the electron layer and the hole layer, 
while the second and third terms describe the exchange contributions from the intralayer interaction and interlayer 
interaction, respectivel y 14 ' 15 . At d <C Ih we have fi^ — \/2n 3 / 2 ^^-n p . Note that the interaction between e — h pairs 
falls off as the distance between the layers decreases, and at d = the e — h pairs form the ideal gas 1 ^. (The last 
statement is valid under the lowest Landau level approximation). 

The nonlinear (cubic in ty) term in the sought-for equation is calculated by substitution of expression ([5]) (without 
derivatives with respect to the "slow" variable R) into cubic in $ terms of the equation Q(ri,r2,t) = 0. Then the 

result should be multiplied by tpo(i) exp^— ^R- (H x f)^ and integrated over f . As a result we arrive at the following 
nonlinear equation for the wave function of the e — h pair condensate: 



ST^ = ~—r SHI^ " — K " |*(IM)|1*(R,t), (13) 



_d V(R.,t) _ h 2 d 2 ^(R,t) //0 
~di ~ ~2M H &R 2 np 

This is nothing but the Gross-Pitaevskii equation in the absence of external fields. From Eq. (fT3|) one can find the 
spectrum of elementary excitations of the e — h pair condensate. It has the well-known form obtained by Bogolyubov 
for a weakly nonideal Bose gas: 




At low k and d ^ the dispersion law is acoustic one, lj = cofc, where the speed of sound is cq = \J /iW /Mh- The 
decrease of the interlayer spacing d leads to the decrease in cq- In the d = case, fi^ = and spectrum (jT4j) become 
quadratic in k. The last, as is known, means the absence of superfluidity. 

From the coefficients entering into Eq. lfTlI)) we can form the parameter £ = h/ ' \/2Mh^ 1 \ which has the dimen- 
sionality of length. This parameter by its sense coincides with the coherence length in the superconductivity theory 
and sets the characteristic scale of wave function module variation. It is not difficult to see that £ ~ l p (lff/d), where 

— 1/2 

l p is the average distance between the pairs (l p ~ n p ). Hence it follows that at d -C Ih the ratio (£/Z p ) 3> 1. In 
this case, a great number of pairs is situated on the scale £, that justifies the macroscopic description of the system 
by the wave function W(R). The coherence length £ becomes of the order of the mean distance between the pairs l p , 
when d ^ Ih, and at d 3> Ih the £ value (formally) becomes much less than l p . This means that in the last case the 
wave function ^(R) loses its macroscopic nature. In other words, there must exist the critical value of the interlayer 
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distance d c ~ Ih , with the excess of which the condition t; l p is violated, and it may be expected that at d > d c 
the Bose condensate collapses. 

Equation (fT"3| can be readily generalized for the case of the electric field E parallel to the layers, provided 
that E weakly varies on the scale Ih- To do this, the potential energies of the electron and the hole (namely: 
Ylkl dr(—l) k eV(r)ip^(r)tpk(^), where V(r) is the scalar potential of the electric field) must be included in the 
Hamiltonian (p]). In this case, the addition eV(r2) — eV(ri) appears in the square brackets of Eq.([7]). Then it is 
necessary to pass on to the coordinates R and f , introduced before, and to exclude the fast variable f by means of 
the above-described procedure. As a result, we come to the following equation 

where a has the meaning of e — h pair polarizability, a(H) — Mh<? / H ■ The equation similar to ()15|) was proposed 
in Ref. 17 on the basis of phenomenological considerations. 



III. VORTICES AND COLLECTIVE EXCITATIONS 



From (fT"5|) it is not difficult to derive the continuity equation d\^\ 2 /dt + divj s = 0, where the superfluid current 
density can be written as j s = I^^Vg. In this case the superfluid velocity v s appears dependent on the electric field 



Here ip is the phase of the wave function of the condensate (\& = | \R | e i</3 ) . 

The right-hand side of Eq.(fT5|) can be considered as a variational derivative 5F/S^>*, where F is the functional of 
the Ginzburg-Landau type. The variation of F in the electric field gives the dipole moment P of the unit area of the 
bilayer system. Using of expression (| 16[) for the superfluid velocity one can obtain the following simple expression for 
P: 



P=«(E+-v s xH)|$| 2 . (17) 

This result means that not only the electric field E, but also the Lorentz force polarizes the pair, acting in opposite 
directions on the positive and negative charges of the pair. 

Eq. (TT5|) has inhomogeneous solutions that describe the quantized vortices. The steady vortex state of the bilayer 
system with vortex centered at points R„ is determined by the relation 



curl *(|j0 =27r^(T n 5(R-R„), (18) 

n 

where a n = ±1 is sign of the n-th vortex circulation, and <5(R) is the two-dimensional S- function. The two-dimensional 
density of the polarization charge (p po i = — divP), associated with the vortices, has the formic 



2-nhH -^-^ 

Ppoi = -an p — — ) er„d(R - R„). (19) 

IVIhC ' — ' 
n 

Hence it follows that vortices have the electric charge q n — ±ev, where v — 2irl 2 H n p . 

In the presence of a uniform electric field the dispersion law of collective excitations of the condensate takes the 
form 



"< k >= k -W(2tt+ 2 <' l,, W < 2 »» 
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where the electric field induced superfluid velocity v s (E) = -jg^-^E x H. Eq. ([20|) represents the dispersion law for the 
excitations propagating in the condensate that moves at a velocity of v s (E). According to the Landau superfluidity 
criterion, the v s value must not exceed the critical value Co = y fi^ /Mh- The critical velocity Co determines the 
critical electric field E cr = (cq/c)H. For n p ~ 10 10 cm -2 , d ~ 10 _6 cm, H ~ IT, e ~ 10, the estimation gives 
E cr ~ 10 3 V/cm. 



IV. BILAYER SYSTEMS WITH CURVED CONDUCTING LAYERS 



So far we have considered idealized structures, where the conducting layers are assumed to be parallel each other. 
In this Section we depart from the idealized situation and take into account a small curvature of the conducting layers. 
In the general case, the equations of curved layers may be written as z — Zi(r), where the index i = 1(2) refers to 
the electron (hole) layer, and r = (x, y) is the two-dimensional radius-vector. The small curvature of the layer means 
a small (in comparison with the interlayer distance) variation of Zi(r) on the magnetic length Ijj scale. Note that 
the curvature of the layers may arise accidentally, but it may also be made purposely, as the structure non-ideality 
leads to a number of rather interesting properties. In particular, the effective mass of the e—h pair, Mh, becomes 
the function of the pair coordinates; the regions arise, where the energy of the e — h pair is lowered (exciton traps), 
and lastly, e — h pairs can gain the polarization in the xy plane. Below we shall demonstrate that this polarization 
can lead to the appearance of quantized vortices even at zero temperature. 

The approach developed in previous sections can be easily extended to the case of curved layers. To do this, the 
argument V12 (i.e., the electron-hole spacing) inEq.([7]) should be written in the form {[Zi(ri) — Z2(r2)] 2 + (ri—r2) 2 } 1 ^ 2 . 
As with the derivation of Eq.QlSp. further transformations come to the transition to the center of mass coordinates of 
the pair R and to the relative coordinate f = ri — r2- Making use of the smoothness of functions Zi(r), we expand 
the expressions Zi(2)(R± r/2) in power f to an accuracy of the second-order terms and then to integrate over the 
fast variable f . As a result, we arrive at the following equation for the wave function of the condensate: 



ih 



d*(R,t)_ W ^d a(R)„ ^ I ( ^ d a(R 

c 
1 



dt 



f^T^r — -E tot x H ) ( ih-^— — -E tot x H 

\ 9R c / 2M H (R) V <9R c 



^a(R)i? 2 ot + C/(R)-M + .9(R)l*(R,t)| 2 *(R,*). (21) 



Here the effective mass Mh of the e — h pair is determined by expression (11), where the interlayer distance depends 
on the coordinate of the pair: d(R) = Z\(R) — ^(R). In this case, the polarizability of the e — h pair a(R) = 
Mh(R.)c 2 / H 2 is also dependent on the coordinates. 

The vector E tot = E + Ej n < where E is the external electric field and E^ is an "internal" field which, just as E 
does, leads to the e — h pair polarization in the xy plane. The field E,- n t is a consequence of the curvature of the 
conducting layers and is given by linear in f terms of the expansion of the functions Zi(R + f/2) and ^(R — f/2): 



e d 
' H 

Here we have introduced the notations Z(R) = (l/2)[Zi(R) + ^(R)] and 



E mt = -^F(d(R)) — Z(R). (22) 




To clarify the physical sense of the field Ei„t, one can imagine a bilayer system, in some region of which the 
conducting layers are the planes parallel to each other, though being deviated from the basal xy plane. In this region, 
the e — h pair has the dipole moment directed along the normal to the conducting layers. Since the normal is deviated 
from the z-axis, the dipole moment acquires the component parallel to the xy plane. The polarization of e — h pairs 
in the xy plane, associated with the conducting layer curvature, can be described through introduction of the internal 
field E int according to Eq. (l22"l) . 

The U(R) term in Eq. (f2Tj) is similar to the potential energy of the boson in the Gross-Pitaevskii equation and is 
determined by the pair binding energy: U(R) = Eo(d(R)). We neglect the corrections to U(R) connecting with the 
finite size of the e — h pair, because they are small as the squared ratio of Ih to the characteristic size of the curvature 
along the xy plane. As concerns the last two terms in Eq. (j21[) . we note that the chemical potential /1 must be derived 
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from the condition of normalization of J |\&| 2 d 2 R to the total number of e — ft, pairs, while the coordinate-dependent 
coupling constant is g(R) — (d(R)) / n p , where the function u^(d) is determined by Eq. (fT2|) . 

It can readily be demonstrated that Eq. (f2"Tj) leads to the continuity equation dn/dt + div(nv s ) = 0, where the 
condensate density n — |$(R, t)\ 2 , and the superfluid velocity is given by: 



fi|l* w xH). (24) 



Mh (R) V dR 

In an inhomogeneous bilayer system, where the effective mass of the pair is dependent on its coordinates, it is 
convenient to pass from the description in terms of superfluid velocity to the description in terms of the generalized 
pair momentum P s = For example, the dynamic equation for the phase of the order parameter, which follows 

directly from Eq. (|2ip , can be written in terms of the generalized momentum of the pair as 



8P S d r 1 ^9 a „ „ ti 2 d / 1 dy / n\ i ,„_. 

^^-M[2lfe P ' + ^ P - E - XH -^M(^^) + ^^ + H- 

One can see from Eq. (|24[) that the difference of the generalized momentum P s from the kinematic momentum p s = 
Mh*v s is connected with the vector E tot x H. An important point is that the circulation of the generalized momentum 
is quantized, while in the general case the circulation of is not quantized. In particular, for the singly connected 
vortex-free region curlP,; = 0, and the kinematic momentum satisfies the equation curlp^ = — S cU v [a(R)E tot ]. In 
other words, even in the absence of quantized vortices, the vector p s may have a vortex component because of the 
inhomogeneous condensate polarization a(R)E tot . Note that in the homogeneous system, where E fot = E, the last 
equality reduces to the known condition curlv s = if divE = 0. Thus the condition curlv s = is a special case of 
the more general quantum condition §V s dR = 2irhk, where k is an integer. 

The dipole moment P of the unit area of the inhomogeneous bilayer system is given by expression (17), where E 
must be replaced by Etot ■ Calculating div P we obtain the following expression for the two-dimensional polarization 
charge density: 



cn c f on \ . . 

p poi = -_curl z P s --(— xP s j z . (26) 

Eq. ({2"6")) can be used to discuss the peculiarities of vortices in inhomogeneous bilayer systems. We restrict the discussion 
to two cases. First we consider the vortex at the center of the axially symmetrical trap for e — h pairs. From Eq. (|26p 
it follows that for this vortex the polarization charge density is given by 



Ppo i = -2irael 2 H n(0) S(R) - ael 2 H - — . (27) 

Here a = ±1 is the sign of the vortex circulation. By integrating Eq. (|27|) over R we obtain that the charge concentrated 
in the region of radius R is determined by the e — h pair density at the boundary of the region: q(R) = —2it<jel 2 H n(R). 
This expression transforms to the formula derived in the previous section for the coordinate- independent vortex charge. 
Indeed, in the homogeneous system, far from the vortex core, the density n(R) quickly approaches the constant n p 
value. Therefore, at R ^> £ the vortex charge is independent of R and is equal to q — —aev. On the other hand, if 
the vortex is located at the center of the trap, where the captured e — h pairs are contained in the R < Rq region, 
then the total electric charge associated with the vortex equals zero (since n(R) = at R > Rq). The contradiction 
with the previous result is removed if one takes into account that in the homogeneous system, along with the vortex 
charge q, there arises the opposite charge (—q) at the boundaries of the system. 

To conclude this Section, we discuss the question whether the curvature of conducting layers may give rise to 
vortices. For simplicity, we consider the axially symmetrical situation, Z — Z(R), and put d(R) — const. In this case, 
we have the effective mass of the e — h pair Mh — const, the field E tot = Ej„t = -^rF(d)^e r , and the superfluid 

velocity v s = ~ ^(^)f^) e 9- Here k is the integer, e r and e# are, respectively, the radial and azimuthal 

unit vectors. The fc-dependent contribution to the energy of condensate is given by the integral J (Mh /2)v 2 \^\ 2 d 2 R. 
At d <C Ih the condensate energy is minimal when 



(28) 
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where Rq is the radius of the system, and the function Int(x) defines the integer closest to x. As follows from Eq.(|28[). 
the fulfillment of a rather rigid condition d\Z(R ) — Z(0)\ > l 2 H log(i? /£) gives k ^ 0, i.e., the curvature of layers may 
give rise to quantized vortices even at zero temperature. 



V. CONCLUDING REMARKS 



In conclusion, we justify, in short, the mean-field approximation used above. As was shown for the first time by 
Shick 18 , in a two-dimensional Bose system in the low-density limit the interaction is renormalized essentially with 
respect to the bare one. In particular, for the two-dimensional bosons modeled as hard dicks of the radius a the 
summation of the infinite series of the ladder diagrams yields the effective interaction in the form Aith 2 / mh\{\ / na 2 ) . 
One can think that for the electron-hole pairs considered a similar logarithmic dependence of the interaction on the 
density of bosons would take place as well. In reality, at d -€i Ih it is not the case. 

The validity of the mean-field approximation at d <SC Ih can be justified as follows. In rarefied systems the effective 
interaction is proportional, in the leading order, to the amplitude of scattering / of two bosons on each other. For the 
bosons with the wavelength of order of the average interparticle distance n -1 / 2 in the case of hard-disk scattering the 
amplitude of scattering is / ~ (lnl/na 2 ) -1 and we arrive at the Shick's result. It is important to emphasize that in 
deriving of that result disks considered as impenetrable ones. But the electron-hole pairs may overlap and it is better 
to consider them as "penetrable" discs. If one models the pairs as discs with the interaction potential equal U = Ua 
inside the discs and {7 = outside the discs one can easily obtain the scattering amplitude 

(29) 

(with the logarithmic accuracy). The first term in the denominator can be evaluated from Eq. (|12j) for the chemical 
potential. At d <C Ih for the radius of the disc a = Ih this term is of order of l 2 H / {And 2 ) . Substituting the density 
in the form v /2irl 2 H and again using Ih as the radius of the disc a we obtain the second term in denominator in the 
form h\2ir/v. In real physical systems v > 1/10 (at smaller v the localization effects become important) and the first 
term is much larger than the second one. Therefore, the effective interaction is proportional Uo and coincides with the 
mean-field result. Thus, at d <C Ih the mean-field approximation is self-consistent and its application is well founded. 
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